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ABSTRACT 

We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 
within the context of the Magnetism in Massive Stars (MiMeS) Large Program at the Canada- 
France-Hawaii Telescope. Thirty-two high signal-to-noise ratio circularly polarised (Stokes 
V) spectra and 13 unpolarised (Stokes /) spectra of HD 148937 were acquired in 2009 and 
2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all 
observations (in both Least-Squares Deconvolved (LSD) mean profiles and individual spectral 
lines). The longitudinal magnetic field inferred from the Stokes V LSD profiles is consistently 
negative, in contrast to the essentially zero field strength measured from the diagnostic null 
profiles. A period search of new and archival equivalent width measurements confirms the 
previously-reported 7.03 d variability period. The variation of equivalent widths is not strictly 
periodic: we present evidence for evolution of the amount or distribution of circumstellar 
plasma. Interpreting the 7.03d period as the stellar rotational period within the context of the 
Oblique Rotator paradigm, we have phased the equivalent widths and longitudinal field mea- 
surements. The longitudinal field measurements show a weak sinusoidal variation of constant 
sign, with extrema out of phase with the Ho' variation by about 0.25 cycles. From our con- 
straint on V sin / < 45 km s"' , we infer that the rotational axis inclination ; < 30°. Modeling the 
longitudinal field phase variation directly, we obtain the magnetic obliquity /? = 38^28° ^^'^ 
dipole polar intensity Z?d = 1020;|^3jq G. Simple modeling of the Ha equivalent width variation 
supports the derived geometry. The inferred magnetic configuration confirms the suggestion 
of Naze et al (2010), who proposed that the weaker variability of HD 148937 as compared 
to other members of this class is a consequence of the stellar geometry. Based on the derived 
magnetic properties and published wind characteristics, we find a wind magnetic confinement 
parameter rjt - 20 and rotation parameter W = 0.12, supporting a picture in which the Ho' 
emission and other line variability have their origin in an oblique, rigidly rotating magneto- 
spheric structure resulting from a magnetically channeled wind. 

Key words: Stars: magnetic fields - Stars : rotation - Stars: early-type - Stars: binaries: 
spectroscopic - Instrumentation : spectropolarimetry. 
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1 INTRODUCTION 

The enigmatic Of?p stars are identified by a number of peculiar and 
outstanding observational properties. The classification was first 
introduced by Walborn (1972) according to the presence of C m 
/i4650 emission with a strength comparable to the neighbouring 
N m lines. Well-studied Of?p stars are now known to exhibit recur- 
rent, and apparently periodic, spectral variations (in B aimer, He i, 
C m and Si ni lines), nan'ow P Cygni or emission components in 
the Balmer lines and He i lines, and UV wind lines weaker than 
those of typical Of supergiants (see Naze et al. 2010 and references 
therein). 

Only 5 Galactic Of?p stars are known (Walborn et al. 

2010) : HD 108, HD 148937, HD 191612, NGC 1624-2 and 
CPD-28°2561. Three of these stars - HD 108, HD 148937 and HD 
191612 - have been studied in detail. In recent years, HD 191612 
and HD 108 have been carefully examined for the presence of mag- 
netic fields (Donati et al. 2006; Martins et al. 2010; Wade et al 

2011) , and both have been clearly detected. In retrospect, the sim- 
ilarity of some of the observational properties of the Of?p stars to 
the 07 V magnetic oblique rotator 6' Ori C appears to have been a 
strong indicator that they are also oblique rotators with strong, sta- 
ble, organised magnetic fields. As the brightest known Of?p star, 
and in light of the report of the marginal detection of a longitudinal 
magnetic field by Hubrig et al. (2008, 201 1), the time is ripe for a 
more intensive investigation of HD 148937. 

HD 148937 has been extensively observed; its observational 
properties are discussed in the literature by Naze et al. (2008a, 
2010). Similar to other well-studied Of?p stars, HD 148937 is re- 
vealed to be a high mass (50-60 M©) main sequence O-type star 
with an effective temperature of about 40 kK. The physical param- 
eters of this star (as first derived by Naze et al. 2008a and revised in 
this paper) are reported in Table 1 . The HR diagram positions of the 
3 well-studied Of?p stars are illustrated in Fig. 1. While it is clear 
that the Of?p stars are high-mass O-type stars, their HR diagram 
positions are relatively uncertain due to their poorly constrained 
distances. This translates into important absolute and relative un- 
certainties on their masses and ages, as is reflected in the relative 
large error bars in Fig. 1 . 

Like HD 108 and HD 191612, HD 148937 is a spectroscopic 
variable star. However, it is distinguished from those stars by its 
lower-amplitude variability, particularly in the primary diagnostic 
Hor line. Whereas HD 108 and HD 191612 exhibit Ho- equivalent 
width (EW) variability of more than a factor of 2, the Ho- line of 
HD 148937 changes in EW by only about 20% (e.g. Naze et al. 
2008a). Despite this much weaker variability, Naze et al. (2008a) 
were able to tentatively identify a period of 7.031 ± 0.003 d in EW 
measurements of HD 148937 spanning over 3 years. This period 
was subsequently confirmed using higher quality data by Naze et al. 
(2010). Photometry of HD 148937 reveals no significant broadband 
flux variation (Balona 1992, van Genderen et al. 1989, 2001, Naze 
et al. 2008a). 

Naze et al. (2008a) also found that this star exhibits an en- 
hanced nitrogen abundance (by about a factor of four relative to 
the sun). At X-ray wavelengths, HD 148937 resembles HD 108 
and HD 191612 in having a thermal spectrum dominated by a rela- 
tively cool component (kT = 0.2 keV), broad lines (> 1700 km s"' ; 
although see detection of narrower features in high-Z lines by Naze 
et al. 2011), and an order-of-magnitude overluminosity compared 
to normal O stars (log[L;;."^''VLBOL] ~ -6). 

In this paper we perform a first detailed investigation of the 
combined magnetic and variability properties of HD 148937 using 
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Figure 1. HR diagram showing the positions of the 3 Of?p stars in which 
magnetic fields have been detected: HD 108, HD 191612 and HD 148937. 
According to its position on this diagram, HD 148937 is the youngest, 
and most massive and luminous of the well-studied Of?p stai's (although 
it should be noted that the detailed positions are subject to significant un- 
certainties in the distances to these stars, as reflected in the error bars). Evo- 
lution tracks and isochrones include rotational mixing, and are from Meynet 

6 Maeder (2005). Adapted from Maitins et al. (2010). 

an extensive high-resolution spectropolarimetric and spectroscopic 
dataset. In Sect. 2 we discuss the data acquired and the methods of 
analysis used. In Sect. 3 we re-examine the physical properties of 
the star, as well as its projected rotational velocity. In Sect. 4 we 
examine the spectral variability, confirming the periodic variability 
of the Ha and other emission and absorption lines, and the period 
derived by Naze et al. (2008a, 2010). In Sect. 5 we analyse in de- 
tail the magnetic data acquired at the Canada-France-Hawaii Tele- 
scope (CFHT). In Sect. 6 we employ the CFHT magnetic data, in 
combination with the stellar physical and rotational parameters, to 
constrain the magnetic field strength and geometry. Finally, in Sect. 

7 we explore the implications of our study, particularly regarding 
the variability and other properties of HD 148937, the confinement 
and structure of its stellar wind, and of the properties of the general 
class of Of?p stars. 



2 OBSERVATIONS 

2.1 ESPaDOnS spectropolarimetry 

Spectropolarimetric observations of HD 148937 were obtained 
using the ESPaDOnS spectropolarimeter at the CFHT in 2009 
and 2010 within the context of the Magnetism in Massive Stars 
(MiMeS) Large Program and Project. Altogether, 32 Stokes V se- 
quences were obtained - two in May 2009, two in September 2009, 
14 during 7 consecutive nights in June 2010, and 14 during 7 con- 
secutive nights in July 2010. In particular, the final observing runs 
were planned so as to sample two complete rotations of the star, 
hypothesising that the 7.03 d period identified by Naze et al. was in 
fact the stellar rotational period. 

Each polarimetric sequence consisted of four individual 
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Table 1. Summary of stellar physical, wind and magnetic properties of 
HD 148937, derived by Naze et al. (2008a) and in this paper The wind 
magnetic confinement parameter 7j, , the rotation parameter W and the char- 
acteristic spin down time r^pi^ are defined and described in Sect. 7. 



Spectral type 
Teff (K) 
log g (cgs) 

R*(Ro) 
vsini (lans"') 
log{L*/Lo) 
M„ (Mo) 


Of?p 

41 000 ± 2000 
4.0 ±0.1 
15.0 ±2.5 
<45 
5.8 ±0.1 
~ 60 


Naze et al. (2008a) 
Naze et al. (2008a) 
Naze et al. (2008a) 
This paper 
This paper 
This paper 


logM/V7(Moyr-l) 
Voo (km s" ' ) 
Fill factor / 


-6 

2600 

0.01 


Naze et al. (2008a) 
Naze et al. (2008a) 
Naze et al. (2008a) 


Bd(G) 
y8(°) 


1020-3^0 
38!^^ 


This paper 
This paper 


1* 
W 


20 

1.2 X 10"' 
1.6 Myr 


This paper 
This paper 
This paper 


EiB - V) 
logiLx /Lb„i) 


0.67 
-6.1 


This paper 

Naze et al. (2008a) 



subexposures taken in different polarimeter configurations. From 
each set of four subexposures we derive a mean Stokes V spectrum 
following the procedure of Donati et al. (1997), ensuring in particu- 
lar that all spurious signatures are removed at first order. Diagnostic 
null polarization spectra (labeled A') are calculated by combining 
the four subexposures in such a way that polarization cancels out, 
allowing us to check that no spurious signals are present in the data 
(see Donati et al. 1997 for more details on how N is defined). All 
frames were processed using the Upena pipeline feeding Libre ES- 
pRIT (Donati et al. 1997), a fully automatic reduction package in- 
stalled at CFHT for optimal extraction of ESPaDOnS spectra. The 
peak signal-to-noise ratios per 1.8 kms"' velocity bin in the re- 
duced spectra range from 60 to nearly 1000, with a median of 690, 
depending on the exposure time and on weather conditions. Due to 
the low declination of HD 148937, the star was generally observed 
at high airmass (2.8 - 3.5) from CFHT. 

The log of CFHT observations is presented in Table 2. 



2.2 FEROS spectroscopy 

In addition to the 32 circular polarisation spectra, 13 high- 
resolution unpolarised spectra were acquired on 10 nights in 2009 
March using the FEROS spectrograph on the MPG/ESO 2.2-m tele- 
scope at La Silla. These spectra cover a useful wavelength range of 
-360-920 nm at a resolving power of R ^ 48, 000, with a S/N of 
~125 at Ho-. The echelle data were extracted and merged using the 
pipeline FEROS Data Reduction System (see Kaufer et al. 1999 for 
details). The log is given in Table 3. 



2.3 Coralie spectroscopy 

Finally, we have re-analysed the 20 echelle spectra of HD 148937 
obtained in 2008 with the Coralie spectrograph on the 1.2m Euler 
Swiss telescope at La Silla (Chile), described by Naze et al. (2010). 



Table 3. Log of Feros spectroscopic measurements of HD 148937. S/N is 
the inferred signal-to-noise ratio per pixel based on scatter on continuum 
pixels near Ha. Phase is computed using the ephemeris expressed by Eq. 
(2). 



HJD 


Phase 


^exp 


S/N 


2454904.8818 


0.966 


60 


138 


2454905.9261 


0.115 


60 


144 


2454907.9267 


0.399 


60 


119 


2454908.9197 


0.541 


60 


124 


2454909.8465 


0.672 


60 


102 


2454910.9096 


0.824 


60 


115 


2454910.9109 


0.824 


60 


125 


2454911.9068 


0.965 


60 


101 


2454912.8990 


0.106 


60 


139 


2454913.8248 


0.238 


60 


118 


2454913.8262 


0.238 


60 


123 


2454914.9080 


0.392 


60 


130 


2454914.9093 


0.392 


60 


131 



3 PHYSICAL AND ROTATIONAL PROPERTIES 

We have re-determined the main stellar properties of HD 148937 
using the atmosphere code CMFGEN (Hillier & Miller 1998). The 
results are basically unchanged compared to the analysis of Naze 
et al. (2008a). The only minor difference is a slightly higher lumi- 
nosity (logL/Lg = 5.8 vs. 5.75). We consider this to be a better 
determination, based on the fit of the UV-optical-IR SED rather 
than just on the V-band magnitude. We also derive the colour ex- 
cess E{B -V) = 0.67 from this fitting process, assuming a classical 
Galactic extension law (Seaton 1979, Howarth 1983). Fits illustrat- 
ing our results are presented in Fig. 2. 

Based on our spectral analysis and comparison with the line 
profiles of HD 108 and HD 191612, we conclude that macrotur- 
bulence is likely the dominant line broadening mechanism. Tests 
performed using He i /14713 and C iv .{.{5801/5812 reveal that 
convolution of our synthetic spectra with a Gaussian profile (mim- 
icking isotropic macroturbulence) gives better results compared to 
pure rotational broadening. In practice, a macroturbulence of about 
50 kms"' gives satisfactory fits. However, many profiles - in par- 
ticular those of He I - have a strong (wind-induced?) asymmetry, 
and those that have only small asymmetries have broadening that 
is dominated by 'turbulence'. The projected rotational velocity is 
an important parameter needed to constrain the stellar and mag- 
netic geometry. To place a stronger constraint on the v sin i of HD 
148937, we focused mainly on the C iv /15801 line as it is reason- 
ably strong, and reasonably symmetric. 

We used a 40 kK, log = 4.0 TLUSTY (Lanz & Hubeny 2003) 
model for the intrinsic line profiles, slightly scaled to match the 
observed line strength. The model profiles were rotationally broad- 
ened with a simple convolution and a linear limb-darkening coef- 
ficient of 0.4 (numerical integration of specific intensities gives re- 
sults that are indistinguishable); isotropic Gaussian turbulence; and 
a Gaussian instrumental response, before remapping to the sam- 
pling rate of the observations, and adding matching noise. Compar- 
ison of models & observations was performed for a range of v sin i 
and turbulence, in both wavelength and Fourier space. 

We find that the broadening is roughly described by the con- 
straint that 

V<r- + (vsini)" 100 km s"', (1) 
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Table 2. Log of CFHT MiMeS observations of HD 148937, including results of the magnetic analysis described in Sect. 3. S/N is peak signal-to-noise ratio 
per 1.8 km s"' pixel in the reduced spectrum (occurring around Htv). Phase is computed using the ephemeris expressed by Eq. (2). ; = Bz/o" is the detection 
significance of the longitudinal magnetic field. 



Stokes V Null N 



CFHT 


HJD 


Exposure 


Peak 


Phase 


<Bz> ± o"B 




<Bz> ± o-R 




Odometer 


-2450000 


time (s) 


S/N 




(G) 




(G) 




1075687 


2454955.04670 


2720 


927 


0.100 


-81 ± 101 


-0.8 


96 ± 101 


1.0 


1076968 


2454959.97316 


2720 


529 


0.800 


-312 ± 187 


-1.7 


-278 ± 187 


-1.5 


1115108 


2455078.72811 


2400 


532 


0.687 


-242 ± 187 


-1.3 


243 ± 187 


1.3 


1115112 


2455078.75798 


2400 


516 


0.692 


-72 ± 205 


-0.4 


180 ±204 


0.9 


1206040 


2455366.87889 


2400 


780 


0.663 


-440 ± 118 


-3.7 


262 ± 116 


2.3 


1206044 


2455366.90878 


2400 


764 


0.667 


-188 ± 126 


-1.5 


242 ± 125 


1.9 


1206203 


2455367.88596 


2400 


764 


0.806 


-247 ± 122 


-2.0 


66 ± 120 


0.6 


1206207 


2455367.91586 


2400 


721 


0.810 


-308 ± 132 


-2.3 


13 ± 132 


0.1 


1206344 


2455368.85844 


2400 


700 


0.944 


-21 ± 139 


-0.2 


5± 137 


0.0 


1206348 


2455368.88837 


2400 


725 


0.949 


-541 ± 135 


-4.0 


-224 ± 135 


-1.7 


1206558 


2455369.88524 


2400 


506 


0.090 


-214 ± 196 


-1.1 


218 ± 198 


1.1 


1206562 


2455369.91519 


2400 


505 


0.095 


-645 ± 200 


-3.2 


-190 ± 199 


-1.0 


1206741 


2455370.84461 


2400 


389 


0.227 


-112 ±245 


-0.5 


347 ± 246 


1.4 


1206745 


2455370.87453 


2400 


421 


0.231 


-229 ± 228 


-1.0 


235 ± 227 


1.0 


1206893 


2455371.84232 


2400 


760 


0.369 


-235 ± 119 


-2.0 


-97 ± 117 


-0.8 


1207032 


2455371.87222 


2400 


765 


0.373 


-293 ± 128 


-2.3 


-22 ± 127 


-0.2 


1207036 


2455372.84803 


2400 


967 


0.512 


-90 ± 124 


-0.7 


53 ± 123 


0.4 


1206897 


2455372.87793 


2400 


782 


0.516 


-201 ± 123 


-1.6 


63 ± 123 


0.5 


1217152 


2455400.77099 


2400 


721 


0.482 


-213 ± 129 


-1.7 


78 ± 129 


0.6 


1217156 


2455400.80176 


2400 


648 


0.487 


-280 ± 145 


-1.9 


258 ± 145 


1.8 


1217445 


2455401.77634 


2400 


58 


0.625 


2246 ±2710 


0.8 


178 ±2739 


0.1 


1217450 


2455401.81689 


2400 


193 


0.631 


-658 ±551 


-1.2 


282 ± 550 


0.5 


1217692 


2455402.77201 


2400 


251 


0.767 


-465 ± 409 


-1.1 


-163 ±409 


-0.4 


1217696 


2455402.80787 


2400 


317 


0.772 


-160 ±319 


-0.5 


576 ±321 


1.8 


1217846 


2455403.77822 


2400 


587 


0.910 


24 ± 168 


0.1 


223 ± 165 


1.4 


1217850 


2455403.80871 


2400 


582 


0.914 


-58 ± 173 


-0.3 


-174 ± 175 


-1.0 


1218069 


2455404.78463 


2400 


677 


0.053 


-105 ± 138 


-0.8 


-328 ± 138 


-2.4 


1218073 


2455404.81527 


2400 


648 


0.057 


-169 ± 139 


-1.2 


-101 ± 138 


-0.7 


1218262 


2455405.76066 


2400 


728 


0.192 


-156 ± 128 


-1.2 


44 ± 127 


0.3 


1218266 


2455405.79129 


2400 


797 


0.196 


-61 ± 113 


-0.5 


140 ± 114 


1.2 


1218487 


2455406.75036 


2400 


744 


0.333 


-76 ± 123 


-0.6 


-92 ± 123 


-0.8 


1218491 


2455406.78074 


2400 


786 


0.337 


20 ± 115 


0.2 


-143 ±113 


-1.3 



where f is the FWHM (2.354(t) of the Gaussian turbulence. Inter- 
estingly, this is essentially identical to the result for HD 191612 
(Howarth et al. 2007). This is remarkable, given the large differ- 
ence in their inferred rotational periods. As suggested by Wade et 
al. (2011), this supports the view that rotational broadening is a 
negligible contributor to the observed line widths. 

The upper limit to vsini from C iv /i5801 is 60 kms"', and 
I'sini = provides an acceptable fit (with ^ = 106 kms"'). This 
is illustrated in Fig. 3. We also explored a variety of weaker lines, 
but even with such high S/N spectra we are limited by S/N issues. 
Nevertheless, in the end we concluded that the N iv /i5226 line 
gives a slightly better upper limit on v sin ; of about 45 km s"' . 

These numbers are compatible with those reported by Naze et 
al. (2008a). The properties of HD 148937 as derived in this section 
are summarised in Table 1 . 



4 SPECTRAL VARIABILITY 

As described by e.g. Naze et al. (2008a), the complex spectrum figure 4. Combined periodogram for all data sets yielding periods near 

of HD 148937 exhibits many similarities to HD 108 and HD 7.03 d. The adopted period is 7.032 ± 0.003 d. 

191612, as well as some important differences. Like HD 108 and 
HD 191612, the C iii /14650 lines are in emission, with strength 
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Figure 2. Left panel: Fit to the optical spectrum #1206203. The black solid line is the observed spectrum, the dashed red line the best fit model (reff = 40000 K, 
logg = 4.0, logL/Lo = 5.8, M = 10"^ Mo/yr, i'„ = 2600 kms"'). The Balmer lines are not reproduced by the TLUSTY models, (see Sect. 3). Right panel: 
Fit of the UV-optical-lR SED. The red solid line is the best fit model (log L/Lq = 5.8 and E{B - V) = 0.67). The black solid line is the FUSE-IUE spectrum. 
The black points correspond to UBVJHK photometry. 
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Figure 3. Fits to the C iv /i5801 line profile for illustrative combinations of rotation and macroturbulence. The left-hand panel shows that a rotationally 
broadened model matching the fine width fails to match the line shape; an additional broadening mechanism is required. The remaining panels show that the 
projected rotational velocity is poorly constrained, because of this significant turbulent broadening. 



comparable to the nearby N lines. However, according to Naze et 
al., unlike those stars, the C in lines are not observed to vary. In fact, 
due to the substantially higher S/N of the new CFHT data, we de- 
tect weak variability in a large number of spectral lines, including 
the C III lines. 

Following Naze et al. (2008a, 2010) we characterise line vari- 
ability using the equivalent width (EW). Before measuring the EW, 
each spectral line was locally re-normalized, and a telluric correc- 
tion algorithm was applied to regions redward of 5790 A. The EWs 
were then obtained by numerically integrating over the line profile. 



The Icr uncertainties were calculated by propagating the individ- 
ual pixel uncertainties in quadrature. Neither the Coralie nor the 
FEROS spectra had pre-computed uncertainties; therefore, a single 
uncertainty value was estimated from the RMS scatter in the con- 
tinuum regions around the line profile and assigned to each pixel. 

The line showing the strongest variability is Ha, with 
a variability amplitude of 20-30%. The variations of other 
photospherically-important lines showing distinct emission contri- 
butions - Hy, He ii ^4686, ^8, He i /15876 - are also statistically 
significant. He i /16678 is an exception - this line is in full emission. 
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Figure 5. Phased EW variations of lines in tlie ESPaDOnS spectra of HD 148937: Balmer lines Ha, Hfi and Hy; He ii /!4686, He i .15876, He i A6618; He i 
/14471, He I /14713; C iv /i5801, C rv /15812; C m /14650 lines; and the DIB located at 5797 A. Symbols/colours distinguish different cycles. Dashed curves 
represent least-squares sinusoidal fits to the combined ESPaDOnS data. 



but varies only weakly (the variation is detected only in the higher- 
S/N ESPaDOnS measurements, and then only marginally). These 
results are consistent with the report of Naze et al. (2010). 

Other lines of He i (e.g. /14471, /14713), while asymmetric, 
generally show no distinct emission. These lines show no signif- 
icant variability. The C iv /15801, /i5812 lines are slightly variable. 
This appears to be the first detection of EW variability of the optical 
C IV lines in an Of?p star. 

We performed period searches using a Lomb-Scargle algo- 
rithm. As a first step, the procedure was applied to each of the 
timeseries of EW measurements of the strongly variable emission 
lines. We then used the method described by Folsom et al. (2008), 
treating each pixel within the line profile as an independent time- 
series of measurements. Each line profile was interpolated onto the 
wavelength grid corresponding to one of the ESPaDOnS spectra. 
Periodograms were constructed for individual pixels in the profile. 
These periodograms were then averaged together, weighted by the 
amplitude of the variability in their respective pixels. The resulting 



weighted mean periodogram for the profile was then examined for 
a best fit period. 

All periodograms of the emission line EWs show signifi- 
cant power near 7d. These periods, which are reproduced in the 
pixel-by-pixel periodograms with somewhat larger uncertainties, 
are all consistent with the best-fit Ua period derived by Naze et 
al. (2008a). To derive a single period characterising the variability 
of HD 148937, the periodograms resulting from the analysis of the 
EWs and line profile variations of the individual lines were all av- 
eraged. The final resulting periodogram, shown in Fig. 4, exhibits 
a single acceptable period at 7.032 ± 0.003 d. This period is con- 
sistent with that reported by Naze et al. (2010). Based on these re- 
sults, for the remainder of the current study we adopt the following 
ephemeris for HD 148937: 

JD«„m.x = 2454588.67(30) + 7.032(3) • E, (2) 

where we have adopted a zero point corresponding to maximum 
Ha emission as derived from a sinusoidal fit to the entire set of 
phased measurements (i.e. ESPaDOnS, FEROS and Coralie). 
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The phased EW measurements are illustrated in Fig. 5 (for il- 
lustrative clarity we show only the higher-quality ESPaDOnS spec- 
tra). The variations of all variable emission lines are in approximate 
phase with the Ha line. Variability with a similar phasing is ob- 
served in the C iii 4650 A lines (which we note are predicted to be 
blended with an O ii line). Interestingly, the EW variations of the 
C IV /i5801 line appears to be shifted by approximately 0.25 cycles 
relative to the variations of the emission lines, although this shift is 
only marginally significant. We also include EWs measured from 
the DIB at 5797 A to illustrate the lack of variability of a reference 
line that is not formed in the environment of the star. 

We find that in many cases the cycle-to-cycle scatter of the 
EWs is larger than the uncertainties attributable to noise, normali- 
sation uncertainties and telluric line contamination. This appears to 
be a consequence of evolution of the shapes of the EW variations 
of some lines. This is especially evident in the Ha and He ii /14686 
lines which exhibit clear changes in amplitude and shape of the EW 
variation between the two well-sampled cycles (illustrated as green 
diamonds and blue circles in Fig. 5). The cycle-to-cycle differences 
in the EW curves of these two lines track each other almost per- 
fectly, and correspond to clear systematic differences in the profile 
shapes (Fig. 6) that are expressed in essentially all lines to some 
degree. As these lines are located throughout the spectrum, these 
differences cannot be a consequence of systematic errors in local 
normalisation. Nor, given the location of some lines (e.g. /t4686) 
well into the blue part of the spectrum, can it be attributed to tel- 
luric lines. We therefore conclude that this variability is intrinsic to 
HD 148937. 

A potential explanation of these differences would be an error 
in the ephemeris expressed by Eq. (2). To check if this could be a 
plausible cause, we varied the period to see if we could improve 
the agreement between the EW variations from different epochs. 
We still find that the 7.032 d period yields the most consistent EW 
variations amongst different epochs and that changing the period 
within the uncertainty had a negligibly small effect on the rela- 
tive phasing of the different epochs. We also compared the pro- 
files from different epochs corresponding to similar EW values but 
with small relative phase offsets to see if the observed profiles are 
a better match than those with different EWs but essentially iden- 
tical phases (e.g. as in Fig. 6). These comparisons show that we 
find larger discrepancies between the profiles with similar EW with 
small relative phase offsets than those shown in Fig. 6, confirming 
that adjusting the ephemeris is not a solution to this phenomenon. 

Based on these results, we conclude that the differences in the 
line profiles are best explained by intrinsic changes in the amount 
or distribution of emitting material with time, i.e. evolution of the 
magnetosphere of HD 148937. It is unclear if this evolution is sec- 
ular or stochastic. This will be discussed further in Sect. 7. 



5 MAGNETIC FIELD 

5.1 Diagnosis from nightly CFHT spectra 

Least-Squares Deconvolution (LSD, Donati et al. 1997) was ap- 
plied to all CFHT observations. In their detection of the magnetic 
field of HD 191612 in 2006, Donati et al. developed and applied 
an LSD line mask containing 12 lines. Given the similarity be- 
tween the spectra of HD 148937 and HD 191612, we began by 
using this line list to extract, for all collected spectra, mean circu- 
lar polarization (LSD Stokes V), mean polarization check (LSD A') 
and mean unpolarized (LSD Stokes /) profiles. All LSD profiles 



were produced on a spectral grid with a velocity bin of 36 km s"', 
providing reasonable sampling of the observed mean profile and 
enhancing the S/N per LSD pixel. The LSD Stokes / profile shows 
a strong extension to high (~ 350 kms"') velocities in the blue 
wing. This asymmetry is similar to that observed in the He i lines 
of HD 148937, and likely is reflective of the inclusion of such lines 
in the mask. An additional dip and contribution to the depression in 
the blue wing is attributed to the DIB located at ~ 579.8 nm, which 
was not included in the line mask but which is blended with C iv 
/i5801 (which is itself included). 

Using the signal detection criteria described by Donati et 
al. (1997), we evaluated the significance of the signal in both the 
Stokes V and A' LSD profiles. We obtain a marginal detection 
(i.e. a false alarm probability (fap) below 10"'') in one spectrum 
(#1206203). In no case is any signal detected in the N profiles ex- 
tracted from the individual spectra. If we coadd the LSD profiles 
of spectra acquired during a single night (typically 2 spectra per 
night), we obtain one definite detection of signal (false alarm prob- 
ability fap < 10"^) in the V profiles (for the coadded profile corre- 
sponding to spectra IDs 1217692 and 1217696), and one marginal 
detection (1218069 and 1218073). An examination of this LSD 
profile shows a weak asymmetry (positive in the red wing, nega- 
tive in the blue wing) that is not present in A'. 

To quantitatively characterise the magnetic field implied by 
our data, from each set of LSD profiles we measured the mean lon- 
gitudinal magnetic field in both V and A' using the first moment 
method (Rees & Semel 1979) as expressed by Eq. (1) of Wade et 
al. (2000), integrating from -280 to -l- 180 km s"' (as did Wade et al. 
2011 forHD 191612). The longitudinal field measured from Stokes 

V is detected significantly (i.e. \z\ = \B(\/<t > 3) in 3 of our obser- 
vations. In no case is the longitudinal field significantly detected in 
A'. The results of this analysis are summarised in Table 1 . 

Due to the relatively large width of the profile, and the lim- 
ited S/N of our observations, our longitudinal field error bars are 
relatively large: from about 100 G, with a median of 138 G. As ex- 
pected, the longitudinal field measured from the A' profiles is con- 
sistent with the absence of any signal in the A' spectrum. The in- 
ferred field values are distributed randomly around zero in a man- 
ner that is statistically consistent with a Gaussian distribution (as 
characterised by the detection significance z). The longitudinal field 
measured from the A' profiles exhibits a mean value of 62 ± 26 G 
(2.4(t). On the other hand, the longitudinal field measured from the 

V profiles is generally negative, with a (weighted) mean value of 
-143 ± 26 G (5.5(t). This again suggests the presence of a signal in 

V that is not detected in A'. 

As a further test, we have investigated if the distribution of 
longitudinal field detection significance from Stokes V is in fact 
significantly different from that derived from N. We employed a 
two-sided Kolmogorov-Smirnov test. For the cumulative distribu- 
tions of zv and zn< we obtain D = 0.5625, indicating that the null 
hypothesis (that the distributions are the same) can be rejected at 
over 99.9% confidence. 

We therefore conclude that the data provide strong evidence of 
the presence of a magnetic field in the photospheric layers of HD 
148937. 

5.2 Diagnosis from coadded CFHT spectra 

Because of the magnitude of the mean longitudinal field (~ 150 G) 
and that it appears to remain consistently negative, it is reasonable 
to assume we are seeing an organised magnetic field, viewed pre- 
dominantly from a single magnetic hemisphere (i.e. near the nega- 
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Figure 6. Cycle-to-cycle changes oftheof Ha, Yip, Hei/i5876, Heii/i4686, He i i4471, C iv ,15801 lines of HD 148937. Shown are profiles obtained at nearly 
identical phases according to Eq. (2). Profiles acquired close in time to each other (#1206203 - phase 0.8080 - and #1206207 - phase 0.8123) match perfectly, 
while a profile obtained 408 days earlier (approximately 58 cycles earlier) at almost identical phase (spectrum #1076968 - phase 0.8008) shows systematic 
differences unattributable to local normalisation. 



tive (south) magnetic pole). In this case, even if the star is rotating 
according to the ~ 7.03 d period, the Stokes V Zeeman signature's 
shape would vary relatively little from one observation to the next. 
We therefore proceed to average all LSD profiles to obtain a single 
grand average profile, with the hope that we are able to detect a 
significant Stokes V signature. 

The grand average LSD profiles are shown in Fig. 7. The A' 
and V profiles are characterised by an LSD noise level (per 36 
kms"' pixel) of 2.4 x 10"^. A clear variation is visible in Stokes 
V, whereas A' shows no variation whatsoever. The detection prob- 
ability of a signal within the line in Stokes V is effectively 100% 
(i.e. fap of 1.1 X 10"^), corresponding to a definite detection. No 
detection is obtained in N (fap of 1.2x 10"'). The longitudinal field 
inferred from the grand average V profile is -193 ± 26 G (7.4(t), 
while no similar field is detected in N (37 ± 26 G; 1 Acr). The detec- 
tion of a clear signature in the grand average profile provides clear 
confirmation that HD 148937 hosts a magnetic field, with a strong 
disc-averaged and phase-averaged longitudinal component. 

We have also searched the grand average spectrum for Zee- 
man signatures in individual line profiles. As illustrated in Fig. 8, 
marginal signatures are visible in Hj8, He i /1541 1 and C iv /t5801, 
while a rather strong signature is visible in He i A5S76. The shapes 
of these signatures are compatible with that of the Stokes V LSD 
profile. Investigation of individual and phase-binned (see below) 
spectra provided no detection of signatures in individual lines - they 
are only detectable in the grand average spectrum. 

To test the sensitivity of these results to the detailed mask com- 
position, we follow the procedure employed by Wade et al. (201 1) 
for HD 191612, extracting LSD profiles for two additional line 
masks. We began using a generic line mask based on a 40000 K ex- 
tract STELLAR request from the VALD database, using a line depth 
threshold of 0. 1 . We sought to construct a mask that yielded a mean 
LSD Stokes / line profile that exhibited as little variability as possi- 
ble, and that was as symmetric as possible, while still maximizing 
Stokes V signal. At each step of the mask development, we visu- 



alised the agreement of the LSD model (i.e. the convolutions of the 
Stokes / and V LSD profiles with the line mask) with the reduced 
spectrum, and evaluated the symmetry and variability of the LSD 
profile. After removal of the most significantly affected lines, we 
were left with a mask containing 26 lines, for which we adjusted 
the mask depths to best match the mean observed line depths. This 
mask is dominated by lines of neutral He, but also contains lines of 
ionised He, C, N O and Si. LSD profiles extracted using this mask 
- which contains more than twice the number of lines in the Donati 
et al. mask for HD 191612 - yielded LSD profiles nearly identi- 
cal to those obtained from the Donati et al. mask. This is likely a 
consequence of the inclusion of many of the same strong He lines 
in both masks. Measurements of the longitudinal fields extracted 
from these LSD profiles are in good agreement with those obtained 
above. 



5.3 Further tests 

As a second test, we extracted LSD profiles using a more restricted 
mask containing only a half-dozen metallic lines (lines of C, O 
and Si). This mask yields profiles whose shapes differ substantially 
from those obtained from the first two masks, as is expected due 
to the exclusion of the much broader He lines. Longitudinal fields 
computed using this mask (measured using an integration range of 
-175 to -1-115 kms"') are in formal agreement with those obtained 
using the masks containing He lines, with a mean of -160 ± 23 G. 
The grand average metallic line LSD profile yields a marginal de- 
tection in Stokes V (and no detection in A'), and corresponds to a 
longitudinal field of -121 ± 23 G (5.3o-). 

As a third test, we compared the nightly-averaged profiles pro- 
duced by the metallic line mask to a grid of synthetic Stokes V 
profiles using the method of Petit & Wade (201 1). We restrict our- 
self to the metallic line mask for this test because in order to treat 
the LSD profiles as a single line profile for modelling, it is best to 
use lines with similar shapes. The synthetic Stokes V profiles were 
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Figure 7. Grand average LSD profiles of HD 148937. From these profiles 
we obtained a definite detection in V, and no detection in N. The longitudi- 
nal field measured from the Stokes V profile within the indicated integration 
range is = -193 ± 26 G (7.4it). For comparison, the longitudinal field 
inferred from the diagnostic N profile is 37 ± 26 G (I Act). The false alarm 
probability in Stokes V is 1.1 X 10"**. while in W it is 1.2 X 10"'. 
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Figure 9. The probability density functions marginalised for have been 
normalized by their peak values in order to facilitate graphical represen- 
tation. The parameter evaluation for the dipole model treat any difference 
with the model as additional Gaussian noise that is marginalized, leading to 
the most conservative estimate of the parameters. 
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Figure 8. Profiles of individual lines in the grand average spectrum of HD 
148937. Note in particular the cleai' Stokes V signature in the absorption 
component of He i i5876 reproduced in two adjacent orders, and a compat- 
ible but much weaker signature in the C iv i5801 line. 

produced by a simple algorithm performing a disk integral and as- 
suming the weak field approximation to compute the local Stokes 

V profiles. The parameters of the synthetic fine were chosen to fit 
the intensity LSD profile. The magnetic model is a simple centered 
dipolar field, parametrized by the dipole field strength B^, the ro- 
tation axis inclination ; with respect to the line of sight, the posi- 
tive magnetic axis obliquity /3 and the rotational phase (p. Assuming 
that only ip may change between different observations of the star, 
the goodness-of-fit of a given rotation-independent (Sj, /3) mag- 
netic configuration can be computed to determine configurations 
that provide good posterior probabilities for all the observed Stokes 

V profiles, in a Bayesian statistic framework. Notice that in order 
to stay general, we do not constrain the phases nor the inclination. 
The Bayesian prior for the inclination is described by a random ori- 



entation (p(i) oc sin(0 di), the prior for the dipolar field strength has 
a modified Jeffreys shape with a cut at = 40 G corresponding to 
twice the grid step in order to avoid a singularity at Bd = G. The 
obliquity and the phases have simple flat priors. In order to assess 
the presence of a dipole-like signal in our observations, we compute 
the odds ratio of the dipole model with the null model (no magnetic 
field, Stokes V = 0). We also perform the same analysis on the null 
profiles. As expected, no single observation is strongly in favour 
of the magnetic model, although most of the observations that do 
favour the magnetic model are clustered around the same phases. 
Taking into account all the observations simultaneously, the odds 
ratio is 22: 1 in favour of the magnetic model. For A', the combined 
odds ratio is 2: 1 in favour of the null model. Note that as the case 
6d = G is included in the magnetic model, the difl'erence between 
the two models for noise is expected to be dominated by the ratio of 
priors, i.e the Occam factor that penalizes the magnetic model for 
its extra complexity. The derived probability distribution functions 
(PDF) marginalised for the dipole strength are shown in Fig. 9. The 
most probable dipole strength according to the metallic line Stokes 
V profiles is about 350 G - approximately 3 times greater than the 
longitudinal field inferred from those profiles. 

Based on this analysis, we conclude that the shapes of the LSD 
profiles depend on the composition of the line mask. However, all 
analyses, independent of the line mask used, imply the existence of 
a magnetic field with a magnitude of the longitudinal component 
of 100-200 G. Given these results, and for consistency with the 
analysis of HD 191612 by Wade et al. (2011), we have based our 
following analysis on the LSD profiles extracted using the 12 line 
mask of Donati et al. (2006). 

5.4 Longitudinal magnetic field variation 

To study the magnetic field geometry, we model the longitudinal 
field variation of HD 148937 as a function of rotational phase. To 
this end, we proceed as follows. First, we adopt the oblique ro- 
tator model as the framework in which we interpret the magnetic 
observations. This implies that the 7.03 d period (and in particular 
the ephemeris as expressed by Eq. 2) corresponds to the stellar ro- 
tational ephemeris. We have binned the LSD profiles in 7 phase 
bins (phases determined according to Eq. 2) typically spanning 
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Table 4. Longitudinal magnetic field of HD 148937. binned in phase ac- 
cording to Eq. (2). P\i and Pn are the detection probabilities derived from 
the LSD Stokes V and N profiles, respectively, z = B^jtr is the detection 
significance of the longitudinal magnetic field. 



Phase 


Stokes V 






NullW 






<Bz> ± o-B 


zv 


Pv(%) 


<Bz> ± o-B 


ZN 












(G) 






0.081 


-189 ±63 


-3.0 


99.906 


-52 ± 63 


-0.8 


04.703 


0.214 


-96 ± 74 


-1.3 


00.129 


156 ±74 


2.1 


82.031 


0.355 


-119±61 


-2.0 


46.803 


-74 ± 60 


-1.2 


96.595 


0.501 


-218 ±65 


-3.3 


70.384 


80 ±65 


1.2 


15.743 


0.662 


-285 ± 74 


-3.8 


99.087 


254 ± 74 


3.4 


96.646 


0.793 


-296 ± 77 


-3.8 


38.587 


8 ±77 


0.1 


66.664 


0.931 


-183 ±75 


-2.4 


99.999 


-50 ± 74 


-0.7 


87.894 




0.03-0.04 cycles (but as large as 0.06 cycles), and re-measured the 
longitudinal field from these averaged LSD profiles. The binned 
V longitudinal field measurements are consistently negative, and 
a majority of them correspond to detections at 3cr confidence or 
greater. In contrast, the associated measurements obtained from the 
N profiles scatter randomly about zero. (While one of the A' lon- 
gitudinal field measurements is formally significant at over 3cr, no 
corresponding coherent signal is apparent in the A' profile). The re- 
sultant uncertainties are reasonably uniform, about 60-75 G. These 
binned results are reported in Table 4. 

The phase variation B({(f>) of the longitudinal field is illustrated 
in Fig. 10. A fit by Least-Squares of a cosine curve of the form 
Bf(0) = Bo + Bi cos(2;r(0 - 0o)) to the data yields a reduced 
of 0.3, with parameters Bq = -200 ± 9 G, B, = 90 ± 12 G and 
00 = 0.24 ± 0. 15. Therefore, according to Least-Squares, the longi- 
tudinal field is significantly different from zero (at over 20cr), and 
the variation of the field is detected at over la. On the other hand, 
the phases of the extrema are relatively weakly constrained. The 
reduced of the data relative to the straight line Bf(0) = (the 
hypothesis of a null field) is 10.2, indicating that the null hypothe- 
sis does not provide an acceptable representation of the data. While 
the reduced x^ of the cosine fit to the data is relatively low, it is not 
unreasonable given the small number of degrees of freedom. 

For illustration, we include the 4 measurements reported by 
Hubrig et al. (2008) and Hubrig et al. (2011) in Fig. 10, for both 
their full spectrum and hydrogen line results. Those measurements 
are qualitatively and statistically consistent with the variation we 
derive. 



Figure 10. Longitudinal magnetic field variation according to the 7.032 d 
period. Square (red) symbols con'espond to the ESPaDOnS observations, 
while other symbols are the measurements reported by Hubiig et al. (2008, 
2011; green triangles = hydrogen lines, blue stars =full spectrum). 
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Figure 11- X landscape of dipole field strength versus obliquity /3 per- 
mitted by the longitudinal field variation of HD 148937, assuming i = 30°. 



5.5 Surface magnetic field geometry 

With the inferred rotational period and derived radius (from Ta- 
ble 1), HD 148937 should have an equatorial rotational velocity 
Ve = 108 kms"'. Based on the upper limit on vsin/ and its un- 
certainty obtained in Sect. 3 (see Table 1), we obtain / < 30°. We 
therefore conclude that HD 148937 is viewed relatively close to a 
rotational pole, a conclusion consistent with the proposal of Naze 
etal. (2010). 

Adopting / = 30° (which ultimately will yield a lower limit on 
the inferred magnetic field strength), we have fit the phase-binned 
longitudinal field measurements with synthetic longitudinal field 
variations corresponding to a grid of oblique dipole surface mag- 
netic fields characterised by the obliquity angle fS and the polar 
field strength Bj. The best-fit model (according to the x~ statistic) 



is characterised by 6^ = - 1020 G and yS = 38° . The Icr uncertainty 
contours (shown in Fig. 11) permit models with Bj in the range 
-700 G to -1400 G and J3 from 20 to 65°. At 3cr the lower limit 
on the inferred dipole strength is still 450 G, i.e. the presence of a 
magnetic field remains confident. We have also performed the same 
procedure using the unbinned longitudinal field measurements, and 
we obtain results that are statistically identical. 

The sum of the inclination and obliquity angles for models 
able to reproduce the observed variation approximately obeys the 
general rule i+ fS ^ 70°. 

The extrema of the best-fit dipole model occur at phases 0.24 
and 0.74, with approximate Icr uncertainties of ±0.15 cycles. On 
the other hand, the emission line variations (in particular those of 
Ha show extrema at phases 0.0 and 0.5. This will be discussed 
further in Sect. 7. 
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6 STELLAR GEOMETRY ACCORDING TO THE Ha 
MODULATION AND HIPPARCOS LIGHTCURVE 

We have also attempted to constrain the stellar geometry using the 
equivalent width variation of the Ho- line. We have used the "toy" 
model of Howarth et al. (2007) which assumes an oblique, rigidly 
rotating, geometrically thin, optically thick Ha disc that introduces 
modulation of the Ho- emission due to its varying projected area as 
the star rotates. The model includes limb darkening and inner/outer 
disc radii that are specified input parameters. 

We have computed a number of grids of models, varying the 
inner and outer disk radii (inner:outer=1.05:1.10, 2.0:2.5, 1.05:2.5 
in units of R,) and linear limb darkening (from 0-1), for rotational 
axis inclination i from 0-90° and disc obliquity angle a (which we 
assume to be equivalent to the magnetic obliquity p in the magnetic 
model) from - 90°, with steps of 1°. 

For each model we evaluated the relative to the observed 
variation, identifying the 100 best models, which define the locus 
of best-fit a as a function of i. Figure 12 illustrates the distribution 
of these best-fit models in the i/a plane. 

We find that the results are very insensitive to disc size, and 
that the modest sensitivity to limb darkening is readily compen- 
sated by selecting different pairs of (a/i) values. The models are 
also highly degenerate with respect to the angles ; and a. This is in 
strong contrast to the results obtained for HD 191612 (Howarth et 
al. 2007, Wade et al. 201 1). Why do these models constrain the ge- 
ometry so poorly compared to 191612? Because the flat-bottomed 
section of 191612's Ha EW variation gives a strong constraint (the 
disk must be seen nearly edge-on for a significant range of phase) - 
a factor which isn't available for 148937. 

Ultimately, we conclude that for ; < 30° (as implied by the 
upper limit on v sin ; and the 7.03 d period) the disc obliquity a > 
10°, and for i = 30°, 10° < o- ^ 25°. These results are compatible 
with the constraint derived from the magnetic data. Although this 
modeling is very approximate, it reinforces the applicability of the 
CRM and supports the interpretation of the emission line variations 
as due to the varying projection angle of a plasma disc confined to 
the magnetic equatorial plane. 

A second potential constraint on the geometry is potentially 
derived from the photometric variability. Naze et al. (2008a) re- 
viewed published analyses of photometric observations of HD 
148937, and performed their own analysis of Hipparcos and Tycho 
photometry of this star. They concluded that no significant change 
or periodicity is present. Here we compare the observed Hipparcos 
Hp photometric variation, converted to Johnson K-band magnitudes 
according to Harmanec (1998), with the predictions of Townsend's 
Monte-Carlo radiative transfer code for simulating light scattering 
in circumstellar envelopes described by Wade et al. (2011). The 
models were developed for HD 191612, computed for wind con- 
finement parameter rj, (defined in the following section) that dif- 
fers somewhat from that of HD 148937 (77, ^ 50 versus 20), but 
with lightcurves computed for the constraint derived from the mag- 
netic field variation of HD 148937, i+/3 -10°. Unfortunately, the 
available photometry is not sufficiently precise to allow us to obtain 
constraints on the geometry. Nevertheless, the predictions demon- 
strate that high-precision photometry and broadband polarimetry 
would enable the determination of the individual values of the an- 
gles I and yS. 



7 DISCUSSION AND CONCLUSIONS 

We have presented new spectroscopic and magnetic measurements 
of the Of?p star HD 148937 based on extensive spectropolari- 
metric monitoring with the ESPaDOnS spectropolarimeter at the 
Canada-France-Hawaii Telescope, and spectroscopic monitoring 
using FEROS at the ESQ 2.2m telescope. The observations and 
their analysis were undertaken within the context of the Magnetism 
in Massive Stars (MiMeS) Project. 

Using the new spectroscopic observations, in tandem with 
archival UV and IR fluxes, we perform a careful re-determination 
of the stellar physical properties, and place a slightly improved up- 
per limit on the projected rotational velocity. 

Analysis of the new and previously published optical spectro- 
scopic observations confirms the variability of a variety of spectral 
lines, and detects variability of others for the first time. Most, and 
probably all, variable lines vary in phase with Ha. (The C iv /15801 
line may appear to vary out of phase by ~ 0.25 cycles; however, 
the formal uncertainly of the phases of the extrema is sufficient 
large that this cannot be robustly established based on our measure- 
ments.) The equivalent widths measured from these lines confirm 
the variation period reported by Naze et al. (2008a, 2010). We find 
that the variability of the line profiles and equivalent widths is not 
strictly periodic, and that the absorption and emission line variabil- 
ity exhibits cycle-to-cycle changes. This is tentatively attributed to 
evolution of the quantity and distribution of emitting material in the 
magnetosphere of HD 148937. Such a phenomenon has also been 
reported for the other well-studied Of?p stars (Howarth et al. 2007, 
Naze et al. 2008b). 

Our new magnetic observations of HD 148937 confirm the ex- 
istence of an organised magnetic field in the photosphere of this 
Of?p star. The measurements are consistent with a longitudinal 
magnetic field of constant sign and weak variability. If we interpret 
the magnetic field and spectral variability within the context of the 
Oblique Rotator Model (CRM), then the variability period (corre- 
sponding to the stellar rotational period) implies an equatorial rota- 
tional velocity of 108 km s"' . To explain the low inferred v sin i, the 
rotational axis of HD 148937 must be viewed at relatively low in- 
clination, i < 30°. Rotational modulation at low inclination is able 
to qualitatively and quantitatively reproduce many of the observed 
properties of HD 148937, in addition to the low vsini: a longitu- 
dinal magnetic field of constant sign and approximately constant 
strength, the single-wave nature, sinusoidal shape and low ampli- 
tude of the Ha EW variation, and the lack of observed photometric 
variability. 

We modelled the longitudinal field variation of HD 148937 
according to the dipole ORM, assuming a rotation axis inclina- 
tion to the line of sight i = 30° (providing a lower limit on the 
strength of the (assumed dipolar) surface magnetic field). We de- 
rive the dipole polar intensity Ba = 1020^3jjJ G and obliquity angle 
= 38^^g°. For other values of the angles follow the approximate 
relation i+f} ^ 70° . Using the stellar and wind parameters reported 
in Table 1 , we compute the wind magnetic confinement parameter 
77, = BI^^R^ / Mv„ 20 (neglecting clumping) and rotation param- 
eter W = Veq/vcrii = 0.12. This places the Alfven radius at about 
^Aif = ^'J* Rt = 2.0 R,. The Kepler (or corotation) radius is located 
further from the star, at about i?Kep = W^^^^ R, = 4.0 R.,. As pointed 
out by ud Doula et al. (2008), for any material trapped on magnetic 
loops inside the Kepler radius, the outward centrifugal support is 
less than the inward pull of gravity; since much of this material 
is compressed into clumps that are too dense to be significantly 
line-driven, it eventually falls back to the star following complex 
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Figure 13. Photometric (left) and linear polaiisation (centre and right) phase variations of HD 148937 computed for vaiious rotational axis incHnations. 
The open circles in the left-hand panel represent the Hipparcos photometry converted to Johnson V band. The Hipparcos photometry of HD 148937 is not 
sufficiently precise to consti'ain the inclination. Nevertheless, it is clear that high-precision flux and polarisation measurements would yield new consti'aints on 
the stellar geometry. 
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Figure 12. Stellar geometry implied by the Ha EW variation, according to 
the model of Howarth et al. (2007). Dift'erent curves coiTespond to dift'ering 
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patterns along the closed field loops. Hence, for HD 148937 (as for 
HD 191612; Wade et al. 2011), all magnetically-confined plasma 
(i.e. all wind plasma located inside the Alfven radius) is unstable to 
this phenomenon. HD 148937 therefore appears to lack a region in 
its circumstellar environment where centrifugal forces are able to 
support magnetically-channelled wind plasma against infall back 
toward the stellar surface (although we recall that the uncertainties 
on T], are quite large, and therefore the location of Rah insecure). 
In any case, while the wind confinement parameter of HD 148937 
is comparable to those of the other known magnetic Of?p stars HD 
191612 and HD 108 and 0' Ori C (~ 50) , its relatively rapid ro- 
tation (P„, ~ 7 d versus ~ 538 d (HD 191612), ~ 55 y (HD 108) 
and 15.5 d (9' Ori C) may lead to a greater influence of rotational 
centrifugal effects on its wind confinement. 

The above discussion implicitly assumes that the variable line 
emission arises in a reasonably static (in the rotating stellar frame) 
distribution of plasma located principally in the magnetic equatorial 



plane. We point out that the He ii A 4686 and H a /i6563 line shapes 
could also be interpreted as P Cyg profiles produced in an outflow. 
If we assume that these profiles arise from a spherically-symmetric 
wind (admittedly a rather rough assumption given the observed line 
variability, and the sensitivity of mass-loss rates determined from 
recombination lines to the (rather uncertain) volume filling factor), 
we can use the CMFGEN code to model these lines using the physi- 
cal parameters in Table 1 . We find that the profiles of these lines can 
be reproduced with an average "density" M/i'oo of about 1.8 x 10"'. 
This is a of factor five larger than that is needed to fit the UV lines 
(M/Voo = 3.8 X 10"'", neglecting clumping). The width of the emis- 
sion line profiles however suggests a substantially lower flow ve- 
locity in the line formation region (~ 900 km s"' ). This results in 
an only moderately higher mass loss rate (1.6 + 0.4 x 10"* Mg/yr) 
than inferred by Naze et al. (1 x 10"* Mo/yr). As Ha is usually 
formed within one stellar radius above the stellar surface (i.e. be- 
low the Alfven radius in the present case), the line emission may 
thus originate from a slower, and thus denser outflow in regions of 
closed magnetic loops. 

One peculiarity of the derived magnetic geometry is its phas- 
ing relative to the emission line variations. In the case of the Of?p 
star HD 191612 (Wade et al. 201 1), extrema of the Ha EW varia- 
tion occur simultaneously with phases of magnetic extrema. In the 
case of HD 148937, however, the magnetic extrema appear to be 
offset from the Ha extrema by about 0.25 cycles. In the context of 
a scattering disc located in the magnetic equatorial plane, this phe- 
nomenon is difficult to explain. However, the offset is only a 1.5-2cr 
effect. Given the low amplitude of the longitudinal field variation, 
we speculate that this is likely an artifact. New, higher-precision 
measurements are required to clarify this issue. 

HD 148937 is distinguished from the other well-studied Of?p 
stars (HD 108 and HD 191612) by its relatively short rotational pe- 
riod. If we assume that angular momentum is lost through the mo- 
ment arm generated by the stellar wind and the magnetic field, the 
characteristic "magnetic braking" spintown time Tspin (Eq. (25) of 
ud Doula et al. 2009) is about 1.6 Myr (assuming the characteristics 
of the UV wind, neglecting clumping), and adopting k = 0.17 (ud 
Doula et al. 2009). The spindown time we compute for HD 148937 
is somewhat longer than those of HD 191612 and HD 108 (respec- 
tively 0.33 Myr and 0.83 Myr using recent physical and magnetic 
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data for those stars). Given the ages of the Of?p stars (2-4 Myr; 
ud Doula et al. 2009), these differences alone are insufficient to ex- 
plain the large range of rotational periods of these stars (assuming 
they all began their lives with similar rotation rates) assuming sim- 
ilar ages. The short rotational period of HD 148937 could therefore 
imply that this star is significantly younger than HD 108 or HD 
191612. 

In closing, we note that HD 148937 is surrounded by the bipo- 
lar nebula NGC6164-5. It is known to be enriched (e.g. Dufour et 
al. 1988) and was therefore proposed to be the result of an giant 
eruption of the star. A larger "bubble" surrounds the system, and 
is attributed to wind expansion in earlier phases of the star's life 
(Fairall et al. 1985). The morphology of NGC6164-5 is axisymmet- 
ric, with two bright "blobs" to the north and south of the star. While 
the origin of the nebula is not well understood, such morphologies 
are often explained by ejection in the same general direction as the 
stellar rotation axis,. However, the magnetic results presented here 
rather suggest that this axis is nearly perpendicular to the plane of 
the sky. A detailed modelling of the nebular ejections should thus 
be undertaken, taking into account the geometry suggested by mag- 
netic measurements. 

The main conclusions of this study are that an organised mag- 
netic field with a dipole surface intensity of approximately 1 kG is 
detected in HD 148937. Within the context of the ORM, the Hor 
period of 7.03 d is interpreted to be the stellar rotational period. 
This is, by more than a factor of two, the shortest rotational period 
found so far for a magnetic O star The stellar radius, rotation pe- 
riod, V sin / upper limit, weak longitudinal magnetic field variation, 
Hq- equivalent width variation and lack of photometric variability 
can all be reconciled with the ORM if the star is viewed close to 
a rotational pole (within 30°). While the detection of a magnetic 
field in HD 148937 is now relatively robust, the parameters of the 
oblique rotator model remain rather uncertain - this is a challenging 
object to observe, and its longitudinal magnetic field intensity and 
variability are both weak. 

Improved geometrical constraints should be obtained by mod- 
eling the Ho- equivalent width variation using more sophisticated 
models, and measuring and modeling the continuum flux and lin- 
ear polarisation variation with high precision. In fact, first measure- 
ments of the Ha line linear polarisation of this star were obtained 
by Vink et al. (2009), but with no detection of any stellar contri- 
bution to the measured polarised flux. (We note however that there 
was no detection either for 6' Ori C, suggesting that higher preci- 
sion measurements would be desirable). 

This is the 3rd Of?p star in which a magnetic field has been 
firmly detected, confirming that the underlying cause of their spec- 
tral peculiarities and variations is almost certainly the magnetic 
field. We therefore propose that the spectroscopically-identified 
Of?p stars are themselves a class of magnetic stars - the first among 
the 0-type stars - whose distinctive spectral peculiarities and vari- 
ability identify them unambiguously as magnetic objects. The test 
of this proposal will be to search for an unambiguously detect mag- 
netic fields in the two newest members of this class, CPD -28 2561 
and NGC 1624-2 (Walborn et al. 2010). 
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